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Within the framework of the NJL model, we investigate the modification of the 
pion damping width in a hot pion gas for temperatures ranging from to 180 
MeV. The pion is found to broaden noticeably at T > 60 MeV. Near the chiral 
phase transition T ~ 180 MeV, the pion width is saturated and amounts to 
70 MeV. The main contribution to the width comes from pion-pion collisions. 
Other contributions are found negligibly small. 
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1. INTRODUCTION 

In the atmosphere of growing interest in the quark-gluon plasma, observations of 
dilepton (electron-positron pair) production in relativistic heavy-ion collisions draw more 
and more attention of physicists making efforts to investigate hadron physics in extreme 
conditions of temperature and density. Experimental data have already been taken, and 
some recent analyses on the dilepton production in Pb+Au (158 GeV/u) collisions are 
available, e. g., from the CERES collaboration PQI2]- 

Insofar as the description of hot and dense media directly from QCD is not yet avail- 
able, the modelling of such processes is an urgent problem. Recently, an attempt was 
made to account for the observed dilepton production rate [T|l2*] . using the simple Bjorken 
scenario of the space-time evolution of an ultrarelativistic heavy-ion collision and the 
vector-dominance model expression for the pion electromagnetic form factor jl]. In jl], it 
was concluded that the experimental dilepton spectrum could not be explained without 
assuming a modification of the p-meson mass and width in medium. 

Besides the p-meson mass and width modification, the same is expected for pions. 
Pions dominate in heavy-ion collisions and it is very important to study how their prop- 
erties change with increasing temperature and density, especially when approaching a 
phase transition in a hot matter. The behavior of the pion in extreme conditions has 
already been investigated in models of the NJL type [S1IS1IZ1E] within the mean-field ap- 
proximation, which is common in NJL [3E],but does not take into account collisions 
in the pion gas. However, as we show in our paper, collisions can give a significant 
contribution to the pion width in extreme conditions. 

For light particles (like pions) the density increases with temperature approximately 
as T 3 . Thus, one expects that the density is large near the supposed phase-transition 
temperature T c , and collisions of particles occur at a much higher rate than in cold mat- 
ter. Collisions lead to shorter lifetimes (or larger widths) of hadronic states at extreme 
environmental conditions. In terms of hadron correlators, an additional imaginary con- 
tribution to the hadron self-energy operator comes from collision integrals [H], thereby 
broadening all particles in the hadron gas. The resulting width is closely related to the 
process of returning a disturbed many-particle system to an equilibrium state: damping. 
Further, the width thus formed is called the damping width. 

There are several approaches to calculate the collision integrals. In our work, we 
follow the prescription given by Kadanoff and Baym [TI]. This way, the contribution 
to the hadron self-energy from collision integrals has a straightforward interpretation: its 
reciprocal is the average time between two consecutive collisions, the lifetime. To estimate 
the average lifetime of a hadron state, one needs to know cross sections for different 
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collision processes averaged over the density of particles with the Bose amplification and 
Pauli suppression taken into account. 

One can try to find the modified self-energy of a particle using a self-consistent func- 
tional formalism, as it is described, e. g., by Hees and Knoll in [T2]. However, in |12j . 
the pion damping width was treated as an external parameter and was not fixed. Our 
purpose is to find an estimate for the pion damping width from a simple quark model and 
investigate possible implications of the pion broadening. We calculate cross sections in the 
framework of the bosonized NJL model with the infrared cut-off jH]. The meson spectral 
functions are chosen to be of the Breit-Wigner form. Then we iterate the equation for 
the width until a self-consistent solution is found. 

The structure of our paper is as follows. In Section 2, we introduce the pion damping 
width. In Section 3, the cross section of the process nn — > mi is derived. The numerical 
results are given in Section 4. The discussion is given in the last section. 



In the vacuum, the lifetime of a particular state is determined by the probability of 
its decay into other states. Thus, a stable particle gains no width in the vacuum at all, 
whereas in dense medium, a particle can strike another one and change its initial state. 
So the lifetime of a particular state is determined by its collision rate. There are also 
inverse processes that can restore the decayed state due to the interaction of particles in 
medium, and then one can find this state again. This process prolongs the lifetime of the 
state. Therefore, the average lifetime of a particular state is determined by the direct and 
inverse processes. 

The in-medium hadron properties are usually investigated in terms of correlators. A 
two-point correlator can then be expressed through a spectral function which is just the 
correlator's imaginary part. For a stable particle, whose correlator has only one pole, the 
spectral function is simply a delta function, while for a state whose energy is spread, the 
spectral function is a continuous function of energy and momentum. 

As discussed in the probability, that, after putting a particle with momentum p 
to the gas at the time t, removing the particle with the same momentum at the time t', 
one can find the gas at the same state as in the beginning, decays as e _r< - i ' - ^, where T is 
a constant determining the decay rate and is usually called the width. If one chooses the 
spectral function in the Breit-Wigner form 



2. PION LIFETIME IN HOT MATTER 




(s - M 2 ) 2 + M 2 r 2 



(1) 
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where M is the mass of the state, the probability indeed decays as e~ r< -''~*- ) (see [TT]). 
Here, the width is a measure of the dispersion of the state energy. In general, the width 
is a function of energy and momentum, and one would finally come to a set of functional 
equations for the width, which are difficult to solve. In our work, we approximate the pion 
width by a constant and obtain a simple equation to be solved by iterations. However, 
the width is supposed to depend on the temperature and density. 

The pion damping width (or lifetime r) can be calculated by the following formula 
from [11]: 

r(p) = r- 1 (p) = S>(p)-S<(p), (2) 



where 




E<(p) = / / {2^)% imA \T\ 2 G>{p l )G<{p,)G<{p A ), (3) 



J pi,p;P3,PA 

Pi Pa Pa 



£>(p) = J J j (27r)% im \T\ 2 G<( Pl )G>(p 3 )G>(p 4 ), (4) 

Pl P3 PA 

in accordance with the notation given in ^Tj. Here, Tis the process amplitude (will be 
given in Section 3), Gf(p) = [1 + n»(p, s», T)]Ai(p 2 ); Gf{p) = m(p, s, T)Ai(p 2 ), with rn 
being the boson occupation numbers 

-l 



7lj(p, s u T) 




exp — — — — 1 



(5) 



and Ai(p 2 ) the spectral function of the ith state. We also use the notations J = J 

Pi 

and Sp uP2 ;p3,p4 = 5(pi +P2—P3 — Pi) in and (jlj), omitting the subscript at P2 to comply 
with the general definition of width (J2J); P2 = P is assumed throughout the rest of the 
paper. The integration is performed in the four- dimensional momentum space over the 
momenta P i, p 3 , p 4 . The indices 1 and 2 correspond to the initial states, while 3 and 4 to 
the final ones. The width T is calculated for pion 2 rested in the heat-bath frame. 
For the inverse lifetime of a pion, one thus obtains 

r- 1 = T = J J dsiV Ie iA K (s 1 ) x 
x[n 7r (pi,Si,T)cr dir *(s;si,s 2 ) - 

-(1 + ^(p!, Slj T))cr inv *(s; si, s 2 )}, (6) 

where v re \ is the relative velocity of particles 1 and 2, and <j dir *(s; Si, s 2 ) and 

are the averaged cross sections for the "direct" and "inverse" processes, respectively, with 

the probability of the final states to be off-mass-sehll taken into account: 

a dh ^*(s; Sl , s 2 ) = J ds 3 J ds A A,{s z )A % {s A )a^\s- s u s 2 , s 3 , s 4 ). (7) 
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(See (jT3j) and (fTHjl for definitions of s and Sj.) 

What one needs then is the cross sections of the processes under investigation. The 
amplitudes of these processes can be calculated in an effective model of pion interaction. 
In the next section, we calculate the amplitudes and cross sections in the framework of 
the NJL model with the infrared cut-off jS]. 



3. PION-PION SCATTERING AMPLITUDE AND CROSS SECTION FROM NJL 

MODEL 



3.1. The Tin — > 7T7T amplitude 



Let us consider scattering of n° on a pion from the medium: 7r°7r° — > — > 7r°7r°, 7r°7r° 



7r + 7r~, ty ty — > tt 7T . We allow also for the lightest scalar isoscalar resonance, the a- 
meson, as an intermediate state (toy — > a — > 7T7t), because of its importance shown in 
various investigations of the pion-pion interaction [31311311111 ■ Here, we use an SU(2) x 
SU(2) chiral quark model of the NJL type P1ITU] where, using the bosonization procedure, 
one obtains the Lagrangian for pions and cr-mesons. The part that contains three- and 
four-particle vertices has the form 



£ int = 2mg a a 3 + 2mg n VZa(27T + 7r + (tt ) 2 ) - g\o 



2-2 2 
7T 



9 * Z -{An + K-K + n- + 47r+7r-(7r ) 2 + (tt ) 4 ) - ^a 4 . 



2 v v , v , / 2 

Here, m is the constituent quark mass (m = gnf^; f n = 93 MeV is the pion weak decay 
constant in the vacuum). The constants g n and g a describe the interaction of the pion 
and a-meson with quarks, respectively. They are related to each other by the equation 



g* = g a VZ. (9) 
The constant Z originates from n-ai transitions and is equal to 

6m2V ' (10) 



Ml 



where M% 1 = 1250 MeV is the mass of ai-meson. The values of the constituent quark 
mass and the constants g n and g a were calculated in jHj. In the vacuum we have: m = 242 
MeV, g~]i = 2.61, and g a = 2.18. Their values at finite temperatures and zero chemical 
potential were calculated in j^j; here we use these results as input data in our calculations. 



6 



In our approach, the total amplitude of consists of a contact term and 

three resonant contributions in the scalar channel (see Figs. 1 a — d) 

T n0n o^ n o n o = -24g 2 n + 3T a (s) + 3%{t) + 3%(u), (11) 

where 



Aglm 2 
M a z -x-iM a T a 

and s, t, u are kinematic invariants: 



w = r ^ , (12) 



s = (P1+P2) 2 = (Ps+Pi) 2 , (13) 
t = (Pi-Ps) 2 = (P2-Pa) 2 , (14) 

U = (pi -p 4 ) 2 = -p 3 ) 2 , (15) 

for which the following identity is satisfied: 

s + t + u = si + s 2 + s 3 + s 4 ; Si=p 2 . (16) 

The momentum p\ corresponds to the impacting pion, while p 2 relates to the pion 
rested in the heat bath. The momenta p% and p^ are those for the particles produced after 
collision. 

For the charged pions in the final state, the amplitudes contain two terms: 

T 7T o n o^ 7r + 7T - = ~8gl + 2T a (s) (17) 

for the process 7r°7r° — > n + n~ (Figs, la and 16); 

^(Vi^Tr^i = —8^ + 2%(t) (18) 

for 7T°7r — > 7r°7r ± (Figs, la and lc). 

In medium, model parameters depend on temperature and density. In our work, as 
only mesons are considered, we restrict ourselves to the case of zero chemical potential 
for which the temperature dependence of the model parameter has been obtained in [H]. 
We use these results in jS] as input data in our calculations. 

3.2. Cross sections 

The differential cross section for a pion-pion scattering is determined by the equation 

da dh |7l 2 (l + n 3 )(l +n A ) 



dt 647rs|pi 



cm 



12 



(19) 
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for the direct and 

dt 647rs|pi >cm | 2 

for the inverse processes. Here, pi jC m is the momentum of pion 1 in the center-of-mass 
frame for pions 1 and 2; n 3 and n 4 are the occupation numbers of produced pions. 
Having integrated over t in the interval defined by the lower and upper limits 

, (a x - s 2 - s 3 + sA 2 2 

Vin — I 2^/s J ~ ^^ 1,Cm ' IP3,cm|J , l/J-J 

, ( Si - S 2 ~ S 3 + SA 2 2 
W — I / ~~ vlPl,cm| — IP3,cm|J , 

we obtain the total cross section. The center-of-mass momenta of pions are determined 
by the equation 

|P«,cm| = \l E lcm - S u i = 1,3, (23) 

where the energies E itCm are defined in the center of mass frame 

s + Si - s 2 

El ' cm - 2V~s ' (24) 

^3,cm = 7T~r ' 

l\J S 

To calculate the occupation numbers in (jTUjl and ()2(Jjl. one needs to know the energies 
of pion 3 and pion 4 in the heat-bath frame 

_ S 2 + S 3 - t _ S 2 + S 4 ~ t 

^ 3 " ^7^' 4 " (26) 

In the case of pion-pion collisions, four processes contribute to the total cross section: 

a(im -> tttt) = cr(7r°7r° -> 7r°7r°) + cr(7r°7r° -> tt+tt") + 

+2a(7r 7r+^^ 7r+). (27) 

The process 7r°7r + — > 7t°tt + occurs at the same rate as 7r°7r~ — > 7r°7r~, so we do not calculate 
them separately; we just put the factor 2 at the last term in (|2~7jl . To obtain the pion 
damping width, we substitute the obtained cross sections into Eqns. (JHJ) and (J7J). 



4. NUMERICAL RESULTS 



Using the definition of the pion lifetime given in Section 2 and the cross sections deter- 
mined in Section 3, we evaluate numerically the pion damping width at the temperatures 
ranging from to 180 MeV. The upper limit on the temperature scale corresponds to the 
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expected transition from the phase with broken chiral symmetry to the symmetric phase. 
The resulting curves are shown on Fig. 2. 

All calculations are performed for a neutral pion rested in the heat bath frame. The 
pion self-energy is approximated by a constant. 

As one can see from Fig. 2, the pion state broadens noticeably in a hot matter, as 
compared to the vacuum state, already at T « 60 MeV. At T — 160, the pion-pion 
scattering accounts for about 80 MeV in the total width; this is the maximum value. 
Then the curve in Fig. 2 turns down; this behavior of the damping width is caused by a 
noticeable decrease of the constant g„ after T = 160 MeV. The cross section is proportional 
to g*, whose value reduces by half at T — 180 MeV, comparing to T — 160 MeV. For 
the temperatures from 160 MeV to 180 MeV, this weakening of the pion-pion interaction 
overpowers the expected increase in collision integrals. 

Besides the pion-pion scattering, we estimated also some other possible contributions 
to the pion width. They come from the following processes: nn — > aa, ira — > na, nn — > a, 
and 7T7T — > qq. All these contributions were calculated separately, i. e., with other modes 
switched off. (This gives the upper limit for a solution of the equation for width, because 
the collision integrals decrease if the spectral functions broaden.) They turned out to be 
small. At T = 180 MeV (which is a little below T c = 186 MeV [8J, ita — ► na gives about 
5 MeV, 7T7T — > aa about 1 MeV, and mi — > a even less than 1 MeV. The decays to quarks 
contribute to the pion width less than 7 MeV. The smallness of these contributions was 
the reason to discard them in our calculations whose purpose is to make a qualitative 
estimate for the pion damping width. 

5. DISCUSSION AND CONCLUSION 

In the framework of the 57/(2) x SU(2) NJL model, the pion damping width was cal- 
culated for the range of temperatures from to 180 MeV. The definition for the damping 
width, given in [11,, was used. A self-consistent method where the pion width does not 
depend on energy and momentum was used to estimate the contributions to the pion 
damping width from pion-pion scattering in a hot and dense matter. Upper limits for 
alternative contributions were estimated qualitatively. It was found that the pion-pion 
scattering is dominant, and other processes that give small contributions can be discarded. 

In our investigation we came to the conclusion that, in a hot gas, the pion spec- 
tral function significantly broadens, while nearing the phase transition point. This may 
have many implications for various processes in a hot and dense matter where pions 
are involved. In particular, this can affect the dilepton production through pion-pion 
annihilation in heavy- ion collisions [3]. 
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Of course, a more systematic and formal approach, like one suggested by Hees and 
Knoll ^2], would be more preferable for a study of processes occurring in a hot matter. 
Nevertheless, our approach, which stems from Quantum Statistical Mechanics, has a clear 
physical interpretation, and it allowed us to make a prediction for the behavior of the pion 
damping width in a hot matter. Furthermore, similar calculations can be done for other 
particles, e. g., for the a- and p-mesons. The a-meson width at those temperatures, at 
which the direct decay a — > nn is suppressed and contributions from collision integrals to 
the a- meson width become noticeable, is of particular interest. Moreover, the diagrams 
with the a-pole play an important role in different processes occurring in a hot matter. 

This work has been supported by RFBR Grant no. 02-02-16194 and Heisenberg- 
Landau program. 
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FIGURE CAPTIONS 



Figure 1. Diagrams contributing to the nn — > Tin amplitude. 
Figure 2. Pion damping width as a function of T. 
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FIGURES 



Figure 1: Diagrams contributing to the nn — > ttty amplitude. 
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Figure 2: Pion damping width as a function of T. 
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